BOUNDARY SHAPE CONTROL OF NAVIER-STOKES EQUATIONS 
AND GEOMETRICAL DESIGN METHOD FOR BLADE'S SURFACE 

IN THE IMPELLER* 

KAITAI Lit, JIAN SU*, AND LIQUAN MEI§ 

Abstract. In this paper A Geometrical Design Method for Blade's surface 9 in the impeller 
is provided here is a solution to a coupling system consisting of the well-known Navier-Stokes 
equations and a four order elliptic boundary value problem . The coupling system is used to describe 
the relations between solutions of Navier-Stokes equations and the geometry of the domain occupied 
by fluids, and also provides new theory and methods for optimal geometric design of the boundary 
of domain mentioned above. This coupling system is the Eular-Lagrange equations of the optimal 
pH ' control problem which is describing a new principle of the geometric design for the blade's surface 

' of an impeller. The control variable is the surface of the blade and the state equations are Navier- 

Stokes equations with mixed boundary conditions in the channel between two blades. The objective 
functional depending on the geometry shape of blade's surface describes the dissipation energy of the 
flow and the power of the impeller. First we prove the existence of a solution of the optimal control 
problem. Then we use a special coordinate system of the Navier-Stokes equations to derive the 
objective functional which depends on the surface explicitly. We also show the weakly continuity 
of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's 
surface. 
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1. Introduction. Blade's shape design for the impeller is driven by the need of 
improving performance and reliability. So far we have not found a geometric design 
' entirely from mathematical point of view. As it is well known that the blade's surface 

^SJ , is a part of the boundary of the flow domain in the impeller. We can use the technique 

CO ' of the boundary geometric control problem for the Navier-Stokes equations to design 

the blade's shape. This idea is motivated by the classical minimal surface which is to 
find a surface spanning on a closed Jordanian curvilinear C such that 

o 

oo . J(3) = Aug inf J(S) 

where J(S') = J J dS is the area of S. 

, In this paper we try to propose a principle for a fully mathematical design of the 

' surface of the blade in an impeller. This principle models a general minimal surface 

by minimizing a functional proposed by us. A key point in this modelling process 
is theoretical rationality and the realization of our design procedure. Using a tensor 
analysis technique we realize this procedure and obtain the Eular-Lagrange equations 
for blade's surface which is a system coupling an elliptic boundary value problem, 
the Navier-Stokes equations and linearized Navier-Stokes equations, and prove the 
existence of solution of the system coupling problem. 

This paper is organized as follows. In section [2] we give the main results of this 
paper. In section [31 we derive the rotating Navier-Stokes equations in the channel 
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in the inipeUer with mixed boundary condition under a new coordinate system. We 
give the minimizing functional problem and derive the Euler-Lagrange equations in 
sectional In section [S] another model to design the blade of the impeller is given. In 
the last section we prove the existence of the solution to the optimal control problem, 
including the existence of the solution of Navier-Stokes Equations with mixed bound- 
ary conditions and the weakly continuous dependence relationships of the solution of 
Navier-Stokes equation with respect to the geometry shape of the blade's surface and 
so on. 

2. Main Results. Suppose (a;\a;^) e D C E'^{2D - EudidianSpace). A 
smooth mapping Q{x^,x^) is the image of a surface. On the other hand, suppose 
that (r, 9, z) is a polar cylindrical coordinate system rotating with impeller's angular 
velocity u . 

{er,e0,k) are the corresponding base vectors, z-axis is the rotating axis of the 
impeller. N is the number of blade and e = n/N. The angle between two successive 
blades is ^ . The flow passage of the impeller is bounded by dfl^ = Tin U To^t U Ff U 
TbU S+ U S-. The middle surface S of the blade is defined as the image of the 
closure of a domain D C where 3? : D ^ is a smooth injective mapping which 
can be expressed by that for any point 3fJ(D) € S 

= x^er + x^Q{x\ x^)e9 + x^kyx = ix\x^) € D, (2.1) 

where Q S C^{D, R) is a smooth function, x — (x^, x'^) is called a Gaussian coordinate 
system on S. It is easy to prove that there exists a family of surfaces with a single 
parameter to cover the domain fl^ defined by the mapping D ^ = {R{x^,x^;£,) : 
y{x\x^) e D}: 

R{x^,x^;0 = x^Cr + x^{e^ + e{x^,x^))ee + x^k, (2.2) 

It is clear that the metric tensor a^^ of S'j is homogenous and nonsingular independent 
of ^, and is given as follows: 

dR 3R 

ac.p = -g^Q^ = Sap + r^Mp: a^det{a^p)^l + r^iel + el)>0, (2.3) 

From this we establish a curvilinear coordinate system (a;^,a;^,^) in , 

{r,0,z)^{x\x^,^):x^=z, x^ = r, ^ ^ e-\e ~ e{x\x^)), (2.4) 
that maps the flow passage domain 

= {R{x\x^,^) = x^er + x^{s^ + Q{x\x^))e0+x^ky{x\x^,S,) e Q} 
into a fixed domain in E^{SD Euclidian Space): 

n = {{x\x^) e D,-l <^ <l} in 
which is independent of Surface 5* of the blade, and Jacobian 



therefore the transformation is nonsingular. 
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Assume that (x^ ,x'^,x^) — (r, 0, z),as well known that corresponding metric 
tensor of is {grr = l,g2'2' = t^,93'3' = = OVi' 7^ j'). According to rule of 

tensor transformation under coordinate transformation we have following calculation 
formulae 

dx^ dx^ 

Substituting (j2.3p into above formulae the metric tensor of in new curvilinear 
coordinate system can be obtain 



gai3 = aaf3, 33/3 = = 6/3, 533 = £ , g = det{gij) = E r . (2.5) 
Through this paper we denote — -§^- Its contravariant components are given by 

= <5"/3^ ^3/3 ^ ^/33 ^ _^-lQ^^ g33^^-2^-2(^^^2|y0|2)^ (2.6) 

where |Vep = 6? + and e„ = 
Fust Model: 

Theorem 2.1. Suppose the Q is a blade's surface defined by (j2.ip . Then G 
is proposed as a solution of following elliptic boundary value problem: 

-g|x(r$^(u;, 6)) + w)r = 0, y{x\x^) e D C 

(2.7) 

combing Navier-Stokes equations and linearized Navier-Stokes equations, where 
{w,p) and {w,p) are solutions of compressible or incompressible rotating Navier- 
Stokes equations p.ip or p.7p and linearized Navier-Stokes equations (3.34) 
respectively and 

K°''^^''{w,e) = 2fir^W°"'6''^^, W"'^ ^ j w^wl^d^, (2.8) 

are defined by (4-11) respectively. 
Variational formulation associated with (|2.7p is given by 

f Find 6 e Vr{D) = {q\q e H^D),q\r„ = Go, |^|ro = 0}, 
such that, Vt? e 

+ 1 lo^^^^i'^' 0)^?^ + r^^{w,w, Q)T]]dx = 0. 
Jit Second Model: 

Theorem 2.2. Suppose the Q is a blade's surface defined by (|2.ip . Then Q 
is proposed as a solution of following elliptic boundary value problem : 



Ko{w)Ae + K^^'iw, e)e,A) + ^^'''^(«;)e,e^ 

+F^{w)ex+ Fo{w,e) ^0, (2-10) 
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where Ko{w), K^^iw), F''''{w), F^{w), Fo{w,Q) are defined by (5.15). The 
variational formulation associated with (5.11) is given by 

' Find e e H}^{D) = {v\v e H\D),v = 6* on 7 dD} 
, such that, V?7 G H^{D) 

I /z3 {[«'o(^«,P, Q)V + "fHw,?, Q)VX - tir^W^^%S^]rujer} dx = 0, 



(2.11) 



where 



(2.12) 



and ^o{w,p,Q),^o{w,p),'^^iw,p),+'^^^{w,p) are defined by (5.10)(5.11). 

3. Rotating Navier-Stokes Equations With Mixed Boundary Condi- 
tions. At first, we consider the three-dimensional rotating Navier-Stokes equations 
in a frame rotating around the axis of a rotating impeller with an angular velocity lo: 



( dp 

— -t- div(pw) = 0, 

pa = Vcr + /, 
dT 

pct,(-^ + uPV jT) — div(KgradT) -l-pdivu; — $ = /i. 



(3.1) 



where p is the density of the fluid, w the velocity of the fluid, h the heat source, 
T the temperature, k the coefficient of heat conductivity, C„ specific heat at constant 
volume, and p, viscosity. Furthermore, the strain rate tensor, stress tensor, dissipative 
function and viscous tensor are given by respectively: 



e^j{w) 

e'^{w) 
a'^ iw,p) 

j\ijkm 



= ^(V^Wj + Vj-Wi); i,j = 1,2,3, 

^ 1 
= 5'^ff^"efo„(u;) = -(V'u;^ + V^u;^), 

= A*^"'=™efe™(w;), $ = A^^'=™ey(u;)e,y(u)), 



(3.2) 



where gij , and g*-' are the covariant and contravariant components of the metric tensor 
of dimensional three Euclidian space in the curvilinear coordinate {x^ , x^, ^) define by 
(|2.4p respectively. 



'^^ ^ ^ dxi dx'' dx' 



(3.3) 



The absolute acceleration of the fluid is given by 



dw 



a = — + {wV)w -|-26i;xw + a;x(a;x R), 



(3.4) 



where lu = cuk is the vector of angular velocity, k the unite vector along axis , and R 
the radium vector of the fluid particle. The flow domain Q,, occupied by the fluids in 
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the channel in the impeller. The boundary dCl^ of flow domain fi^ consists of inflow 
boundary F^n, out flow boundary Tout, positive blade's surface 5+, negative blade's 
surface S- and top wall Tt and Bottom wall Ff,: 

dn, = T = F,„ U Tout US-US+UTtUTb (3.5) 

Boundary conditions are given by 

Hs_us+=0, wIrt = 0, w|rt = 0, 

cr*-'(w;,p)nj|r,„ = gin, ^r''' ("^1?)"-^ |r„„t = ^^^^ (Natural conditions) (3.6) 
1^ + A(r - To) = Owhere A > is constant, 

If the fluid is incompressible and flow is stationary then 

diYw = 0, 

{wV)w + 2uj X w + Vp — j/div(e(w)) = —uj x {uj x R) + f, 
w|ro=0, Fo = 5*+ U S*- U Ff U Ffc, 
{-pn + 2i^e(ti;))|r.„ = gin, Ti = Tin U Tout, 
{-pn + 2;/e(i«))|r„„t = gout. 

For the polytropic ideal gas and flow is stationary, system p.ip turns to the 
conservation form 

' div(pu;) = 0, 
div(pu; (g) w) + 2puj xw + BX/{pT) 

= ^Aw + (A + /i)Vdivw - puj X {uj X R), (3.8) 
div[p(J^ + c„T + RT)w] 

= kAT + Adiv(wdivw) + ^div[wVw] + f A|wp, 

while for isentropic ideas gases, it turns 



div(pK;) = 0, 

div(p?« (g) w) + 2puj X w + aV(p^) 

= 2^div(e) + AVdiv7« — pu x {lj x R), 



(3.9) 



where 7 > lis the specific heat radio and a a positive constant. 

The rate of work done by the impeller and global dissipative energy are given by 



I{S,wiS))= / / cr -n-eeUjrdS, J{S,wiS)) ^ / / / ^{w)dV (3.10) 
J Js-us+ J J Jn, 

where eg is base vector along the angular direction in a polar cylindrical coordinate 
system. 

Let us employ new coordinate system defined by ()2.2p . Flow's domain fJ^ is 
mapped into fl = D x [—1, 1], where D is a domain in {x^, x'^) G 3?^ surround by four 
arc AB, CD, CB, DA such that 

9L»==7oU7i, 7o = ABU(7B, -ii^CBiJDA, 
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and there exist four positive functions 70(2:), 70(2;), 71(2:), 71 (z) such that 

r := = 7o(x^) = 7o(z) on AB, =jo{x^) onCD 

r = 7i(x^) = 7i(z) onDA, x"^ —ji{x^) onBC, 

ro < 7o(z) < ri on tq < 70(2;) < ri onCD, 

^0 ^ 7i(^) !i onDA^ rg < 71(2:) < Ti onBC. 



Let 



Ti = Tout U r™, To = U Ft U {e = 1} U {C = -1}, 



(3.11) 



(3.12) 



ai? = 7oU7i, 9r!=jouri, 

70 - (i? n Tfe) u (i? n ro, 71 - u Tout) u (z? u r,„o, 



(3.13) 



where 3? is defined by p.l 
Let denote 



V{n):^{v\,v e H\nf,v\f^^O},H^{n) = {q\,q e H\n),q\f^^Q}. (3.14) 

The variational formulations for Navier-Stokes problem (|3.7p and p.9p are respectively 
given by 



Find(w,p),w G y(f^),pe L2(f]), such that 
a{w, v) + 2{u! X w,v) + b(w, w, v) + 

-(p,divw) = {F,v), yv e V{Q), 
{q,divw) = 0, Vq e L^{n), 



(3.15) 



and 



Find(ui,p),w € V{n),pe L''(f}), such that 
a(w, v) + 2{u! X w,v) + b{pw, w, v) + 

+ (-ap + Adivw;,divu) = {F,v), Vw G T^(f^), 
iVq,pw) = {pwn,q)\r,, Vg e H^i^l), 



(3.16) 



where 



:=< f,v> + < g,v>^^, 
a{w,v) = A'^'''^eij{w)ekTn{v)y/gdxd£,, 



'y/gclxd^, 



(3.17) 



Next we rewrite (jSTT)) and (|3.9p in new coordinate system. Because second kind of 
Christoffel symbols in new coordinate system are 



33 



(3.18) 



p3 
^ 33 
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the covariant derivatives of the velocity field ViW^ = + T^/^w'^ can be expressed 
as 

Lemma 3.1. Under the curvilinear coordinate system (x^,a;^,^) defined by ()2.4p . 
the covariant derivatives of the velocity field can be expressed as 

' Vaw^ ^^-r6!^e^u{w,e), 

VoU;3 = 1^ + s~^{x^)~^w^ea + e-^w/^Oaf} + {ex^)-^a2o.Il{w, 6), 
S/sW' ^^-x^eS2c.n{w,e), (3.19) 

^ divw = ^ + ^ + ^, U{w,e)=ew^ + w^Q(,. 
and the deformation tensors are given by 

^,,iw, e) = v.^(«;)eA + i^^f{w)exe, + e*^^, e), 



where 



r 



, dur 



and 



where 



A^^'''(pki{w)(pij{v)red^dx = {{w,v)) + / B{w,v,Q)d^dx, 

n 



(3.20) 



^3aH =-(_+eV2 — ), ^33(«;)=eV2(— + — ), 



(3.22) 



e;^(u;, 9) = ^r^w^ie^e^p + epQ^a.), el^{w) = i^^w'^e^^, elsiw) = 0. (3.23) 

The proof is omitted here. 

By simply tensor calculations in terms of (|3.19[) (|2.5|) and (|2.6[) . 

A^J'='efci(t/;)e,,(«) 

= e«/3Mea/3(w) + 2(g335"/3 + e-^eaOp)e3aiw)e3p{v) 
+g^^g^^e33{w)e33{v) - 2e^^Q0{ea/3{w)e3a{v) + eap{v)e3a{w)) 
-2e-^g^^ea{e33{w)e3a{v) + e33{v)e3a{w)) 

+e~'^QaQi3{eaf3{w)e33{v) + ecf3{v)e33{w)) = A'-^^'-Lpki{w)^pij{v) 

+A'3^^[ipkl{w)i^^j{v, e) + Vfei(w, Q)f^]{v) + i^kl{w, Q)'l|J^j{v, 6)], 



(3.24) 



{{w,v)) = / 2^ ['Pai3{w)(pap{v) + 2(re) ^<P3a(w)(p3;3(w) + [re) ^'P3'i{w)'P3.3{v)\ red^dx, 
n 

\\w\\h ■■= iiw,w)), 

a{w, v) = {{w, v)) + J2v [B{w, v, 8) + A{w, v, &)] red^dx, 
n 

(3.25) 
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and 

f B{w,v,e) =2Ai[2e-2(|ve|2,5,/3 + e„e0)^3aH^3/3M 

+e-^r-^\Ve\^il + a)(^33Hv33(w) 

-2e~'^g^^ea{(p33{w)ip3a{v) + (p33{v)'P3a{w)) 
+e^'^Qa'S>f3i'Paf3iw)ip33{v) + (fi a (3 {v)(p 33 (w))], 

A{w, v,e): = A'^'''[ipki{w)^,j{v, e) + ip,j{v)^ki{w, 9) + ^J,j{w, e)4>tjiv, 9)] 

(3.26) 

In particular, 

f Biw,w,Q) ^2^[2e-\\Ve\^Sa,p + ea,efj)ip3c.{w)ip3p{w) 

+e-V-2|V9|2(l + a)ip33{w)ip33{w) - 4e-'^Qpipa,p{w)ip3a,{w) 
-4:e~^g^^<da<P33{'w)(p3a{w) + 2e~'^<da'dpfap{w)(p33{w)], 
_ Aiw, w, 9) = A'^'''[2ipklH^^J{w, 9) + ^Pijiw, 9)^„(w;, 9)] 



Next we consider trilinear form. Using (j3.19|) . we have 

w^VjUjI^ = u;"|f^ +^3^ -r(52/3n(w,9)n(w,9), 

+£-V92n(w'i'9)n(u;, 9) + e-^r-^w^iw^ea + n(w, 9)), 



(3.27) 



Therefore 

6(w, w, v) 

= Id I-i 9knrW^^jwH^^gd£,dx = j\{[a^p{yj^^ + w^^) 
+er^Qp{w^^ + + ^)) - r52pT\{w, 9)n(u;, 9) 

+r'^w''w^QaxQ(j + r9^w2(w^9A + n(w, Q))]vf^ 
+er^[Q^{w^^ + w^^) + e{w^^ + + h!)) 

+w°'w'^Qafi + w2(u;"9„ + n(?z;, &))]v'^}redidx 

By similar manner, angular velocity vector is given by 

w = Loei — £^"'^ti'9ie3, 
2uj X w 

= -2ruj5a2{ew^ + w^Qi3)ea + 2rwe~'^[<c)2{ew^ + w^Qp) 
= -2ruj5a2R{w, 9)e^ + 2rwe-i[92n(w, 9) + r-'^w'^]e3, 

while Coliali form 



(3.29) 



C(u;,w) 

— Id /-I ^ wYv^^d^dx 

= Jjj /^J2aa/3(w X + 2e{x'^)^Qa{{c3 x 'u;)"i;3 + (w x wfv"") (3.30) 

+2e^(x^)^(w X w)3w3]ea;^d^da; 

= 2rw[(w29/3 - S2fslliw, e))v'^ + ew^v^]red^dx 

It is very easy to verify that 

C{w,w)=0. (3.31) 
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Throughout this paper, Latin indices and exponents i,j,k--- vary in the set 
{1,2,3}, while Greek indices and exponents a,/?. 7--- vary in the set {1,2} . Fur- 
thermore, the summation convention with respect to repeated indices or exponents is 
systematically used in conjunction with this rule. 

Lemma 3.2. The function \\ ■ \\q defined by (3.25) is a norm in Hilbert space 

V{n):={vG H\Df,v\f^=0,} (3.32) 

Proof. Indeed it is enough to prove that 

= 0,Vu; e w = 0. 

This means that 

\\w\\n = 0,i.e., (pij{w) = 0. 
we have to prove w = 0. Firstly, the following identity is held 

This shows that 

ip„i3{w) = 0, in£) => d-ydaw'^ = 0, mV'{D). 

By a classical result from distribution theory, each function w is therefore a polynomial 
of degree< 1 (recall that the set D is connect). In other words, there exist constants 
Ca and da^ such that 

{x) = Ca + daffx'^ , y X = {x^ , x"^) G D, 

But ^afii'w) = also implies that rfa/j = —dpa\ hence there exist two vectors c,d& 
such that 

w = c + dx x,yx e D, 

Since w;|p = and the set where such a vector field w" vanishes is always of zero area 
unless c = d = 0, it follows that w" = when area Fq > 0. On the other hand, in 
view of boundary condition (3.13) 

at, r at. 

The proof is complete. □ 

Lemma 3.3. The norm \\ ■ ||n and the norm 

i=l a=l ' 

are equivalent in V{0),i.e. there exist a constant Ci{Q) > 0, i = 1,2 depending upon 
only such that 

Ci(0)|«;|i,a < ||«;||q < C2(^^)k|i,n, e V{^), (3.33) 
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Proof. Firstly wc indicate that in view of (3.11)(3.12) we assert that there exist 
a constant Ci{Q) > 0, i = 1, 2 depending upon Q only such that 

3 3 
i,j=l i,j=i 

and ipij{w) can be looked as strain tensor in Cartesian coordinates in 9?^ then ac- 

3 

cording to Korn's inequality (see [14][15]) the ( ^ ll9'ii(^)llo o)^^^ equivalent to 

therefor this reach to (3.33). The proof is complete. □ 
Theorem 3.1. Under the new coordinate system, the stationary Navier-Stokes 
equations can be explicitly expressed as 6; 



dzv{w)^^ + % + ^=0, 

J\f'{w,p,@) ■.= £'{w,p,e)+AfAf\u,u)= f, i=l,2,3, ^" ' 



where 



C"{w,p,e):=-~iy[Aw"-2e-^&0§^ + (re) "a^ + r 

-^rsS^aOx^ - (re)-HSaxe2 + ^S^c&x + rj^x^Q)^ - 2r-M2„^ 
-2e52aaQ2w'' + 52a{r"^52a - 2626^ - a2a|Ve|2 - rQx@Xa)w^)] 
-2r62oc'^Tl{w, 6) + VaP - e-^e^lf, 

AfAf^iw, w) —w^^+w^^- rS2aTi{w, Q)U{w, 6), 



(3.35) 



CHw,p, e) := -,.[Aw^ - 2e-^Qp^ + {re)-^a^ + 3r-'S20^ 
+e~^r-^a^ +2{re)-\S2Bex + rQpx)^ 

-{re)-\Q2 + rAG)^ + (re)-i(r-M2„e2 + 3e2a + rd„Ae)w^] ,3 
-2ra;£-in(«;, 6)62 - 2(er)-ii«2 - £-16/35/3}?+ (re)-ia§|, 
MAf^iw, w) := wl^^ + + e-'^wl^w^e/3x + [rsy^Uiw, e){2w'^ 

+r^Q2liiw,e)), " 

where A(p = -^r§^+-^T§^ , s"^ = +1, —1, depending on {a, /3) = (1, 2), (2, 1), or other case. 
Under the new coordinate system the domain ^ = {{x^,x'^) € D, — 1 < ^ < 1} «s in- 
dependent of geometry of blade's surface G. 

Proof. Substituting (3.19)(2.5)(2.6) and (3.20) into (3.7), tensor calculations show 
that (3.14) is valid. The details is omitted here. □ 

By using above formulae we claim that there exist of the Gateaux derivative of 
the solution of Navier-Stokes equations to be satisfy following linearized Navier-Stokes 
equation 

Theorem 3.2. Assume that there exists a solution (w(G),p(G)) of Navier-Stokes 
problem (3.7) such that define a mapping G (w(G),p(G)) from Hq{D) n H'^{D) 
to H^''^{n) X L'^''^{Q). Then there exist the Gateaux derivatives of {w,p) at a point 

G e H^{D)r\H^{D) with respect to any direction rj € H^{D)r\H'^{D):w = ^^V,P = 
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Vp 

"^^qV '^''^d satisfy following linearized equations: 



div{w) = 

a{w,p,e)+AfAf\w,w) +N'N'\w,w) = R'{w,p) (3.37) 
■f^'lr, =0, a'^{w,p)nj\ro =0- 



where 



+S2^w-i2S2fie, + 2e2h<r + 2r2(e2<5,/3 + '52^©.)|ve|2 

-2a2aQ0 - rQ0a)] - ^[2r52a{Hw, e) + i^)w^ + e-Ho,fi^\ 

(3.38) 

R\w,p, 9) = ve-\-^{^) + id2^{r-hv^)) + e-^dpx{w^w^) - 2u^{Qfi^) 
-u{e)-^^[2{rer^<dp^ + 2r-i^ - r-^^s/S^ + r'^JaX] 
+e~^^[2wP^-2Loruj52pIi{w,@) 

+r{n{w, 6) - uj){2@2W0 - 52pn{w, 6))] + e-^Ap + 2e-^-^{r@0%) 

(3.39) 

Proof The Navier-Stokes equations (3.24) read 



Set Gateaux derivative with respect with 6 along any director G W := H'^{D) fl 
Hq{D) denoted by Then from (3.40) we obtain 

^j\f" {w,p, e)e„77 + (w,, p, 6)637? +^^"iw,p,e)^v+^fHw,p,e)^r] 

^U"{w,p,Q)e^ + ^Af^w,p,Q)es + [Af-{w,p,e) - + [AfHw,P,Q) - /'lifes = 0, 



Hence, 

^U°'{w,p,e)T] = C"{w,%Q)T]+AfAf''{w,w)'n+J\fJ\f°'{w,w)T] + R"{w,p,e)r] = 0, 
^Af^{w,p, e)r] = £^{w,p, &)t] + +AfAf^{w,w)r]+AfAf^{w, w)r] + R^{w,p, 9)// = 0, 



where R'^{iu,p, @) can be obtain from (3.22) (3.23). Indeed, direction calculations from 
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(3.22) show 

+(52a(-2r?2e, - 292//, - 2r\e2Va + %e,)|Ve|2 - 2a2ae0Vi3 - rrixQxa - rQxrixa)w^)] 
-[r52a{mw, 9) + 2uj)wP + w^niw, 9)) + e-H^0^]i^0 

= u{-r62aW''exVXa+e-'^An) + y[2e-^^ - 2{re)-\^@0^ + 2re52c.^]m 
+i/[(r£)-i(5„A^2/3 + 252<^5x0)^ + 2e52ocw''{a520jV 2r''@0@2)]m 

+y52a{'252f3Qa + 292(5/3. + 2r2 (92(5,9 + <52/39,)| VO^ - 2a2ae^ - r@f)„)w'')r]0 

-[2r52cc{Tl{w,Q)+w)wP +e-H^p%]rip 



Taking into accounts of the homogenous boundary conditions for the r} and inte- 
grating by parts and using Green formula and incompressible condition 

1 ~ dw"' „ 1 , ~ dw" 



dx0dx0d£, d^ 



we claim 



i?"(z«,p, 9) = ,yS2ao^{rw''&x) + ^^'^ + ^^Wt'^ ' '^^'^^0^ + '^re52a^' 
+y^[{re)-^So,xS20 + 262o.5x0)^ + 2s52o.w^ {a520 + 2r2e^e2)] 
+J^(52a9|T [(2(52/39, + 292(5/3, + 2r^{Q2Sa0 + (52/39,)|V9p - 2a2,9/3 - r9^,)w;^] 
-5f^[2r<52«(n(w;, 9) + uj)wP + e-H^0%] 



This is (3.25). Next we consider (3.26). Indeed, 

-RS = -^[-2e-iry/3g^ + {rs)-^2r^@f,rjp^ + e-\-^2r^e0n0^ 
+2{re)-H62f3Vx+rvpx)^ 

-{re)-\r]2+rAr])^ + (re)-'^ {r-'^S2aV2 + 5v2a + rdaATj)w°'] + e-'^w^w^rj/sx 
+(re)-iw^?7/3(2w2 + r'^Q^Uiw, 9) - 2w(a;2)292) 

+{re)-^n{w, e){r^r]2ll{w, 9) + r^e2W^rjp - 2uj{x''fi^2) - s-^v^d^p + {re)-^2r'^ep-np^ 
= -zy(e)-i[2|^r7,,A - ^Ary + 'ir-^w"r,2c. + w^d^Ar^] + e-^wl^w\px 
+ye-'[2§^ - 28-^0/,^ - 2{rs)-'Q0^ - 2r-'^ + r-^52p^ - r-H2pw^]v0 
+(re)-i [{2w^ + r'^Q2^{w, 9) - 2ujr'^Q2)wl^ 

+Il{w,ve){r''52i3Tl{w,Q)+r'^Q2wP - 2ijrH20) - rd/3P+2r^e/3^]r]fs, 
By similar manner and 

2 dw^ ~ dw^ 

'^'^W = ^d^ 
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we assert 

+£-1 J^[r«;/' (2^ +262(0(^,6) -w)) +r52/3n(«;,e)(n(«;, 6) - 2a;] 
+e-'Ap + 2e-^^{rQp%) 

= ,ys-\-A(:^) + 3d2a.{r-'w''))+e-'d0x{wf'w^) - 2v^{Qp^) 
-,[e)'^^[2{re)-^Qp^ + 2r-^^ - + r~^S,,w^] 

+£-1 J^[™^(2^ + 2e2(n(u;, e) - w)) + rS20U{w, e){U{w, Q) - 2iv] 
+£-iAp + 2£-igf^(re^||) 



This is (3.26). The proof is complete. □ 

Associational variational formulation is given by 

' Find {w,p) e H^{Q) x ^^(n) such that 
a(w;, v) + b{w; w, v) + b(w; w, v) + 2{lu x w, v) — (p, div?;) , , 

= J^R{w,p)vdV := {R{w,p),v), W v € H^{n), ^ ' 

(g, divw) = o, Vg e L'^{Ct), 

where 

{R{w,p), v) = gijR'v^dV = ^[a„^i?"(«;,p, e)^;'' 

+sr^{R°'{w,p, Q)v^ + R^{w,p, @)v") + e^r^R^{w,p, e)v^]redxd^, 

(3.42) 

For the compressible case, 

div(wp + wp) = 

div(/9w;*w + pw^w + pww^) + 2p{Lu x wY + 2p{ujL0 x w;)* + ag^^V j{'yp'^~^p) 
-S/j{A'^''^eum{w)) = S'{w,p) 

(3.43) 

where 

S\w,p,Q) = -d0S-'P{w,p;@) + dl,S-^^^^-\w,p;Q), ^^^^^^ 



where 

5«''5(w,p, 6) = {r5^[{5{@^ + 5^ex)w^w'' + 2ew^w^] + 2pru>6„2W^ - e-^^a/j^^ 

+2i.e-i[-e„,g^'"V,eL(«') - (e-^QaQa + S,.ijg^^)V 3e:,sH 

-Sap'^xeaxiw) - Vf3eaa(w) - VsCapiw) 
-r~^{Sa2Spx + iSaxSp2){e3x{w) - e"^OAe33(w))]} 
5"'(^^-)(«;,p, ; 6) = -2iye-'[{5a(i6^x + Sa^S0x){e3^M - e-'Q^ea^iw)) 
+e«5^™(V,e^^(«;,e)] 



(3.45) 
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+2r((5f 62 + S!^ex)w^w^ + red!^w^w^] 

+2nj{{w^ + e-^w^ex)Sp2 + e-^w'^e2) - e'^V pp + 2e-'^Qp^ 

+e-^{250x^^e3a{w) - Vxespiw) + V0e3x{w))Q\ 

-2e-2(2V3e3A(«^) + VAeaaH)] - 2i^[-{re)-^e2p{w) 
+e3^{w)[e-^r-\5^2ef3 - ^6020^) + re-\5p^e2 - <52/3e^)|Ve|2] 
+e3s{w){re)-H20{4 + 3r^\VQm 

; e) - {2iy[e-\e-^exea + 5xa9^^)e33{w) - e-^QxSj, + Q^Sxa)e3-yiw)] 
- + 2iyg^^^g^"^Vjei'^{w,e)}. 

(3.46) 



-2„„A„ 



4. A Principle of Geometric Design for Blade's Surface. We consider 

naturcdly to choice global dissipativc energy to be object ftmctional for geometric 
design of the blade surface. The global dissipative energy functional is given by 

J{S) = 1/2/ / / <f{wiS),w{S))dV = 11 J A'^'^'eu{w)e.j{wy9dxd^, 

j^ijki ^ Xg'^g^^ + ^l{g'^g^^ + g'^g'-i^), 

(4.1) 

where D x [— 1, 1] and f2e is flow passage in Turbomachinery bounded by Ft U 
Fft U Tin U Tout U 5+ U S"- . We proposal a principle for geometric design of the blade: 

{Find a surface 3 of the blade such that 
J(9)= m^J(^), (4-2) 

where !F denotes a set of regular surfaces spanning on a given Jordannian curve 
C G . The S is called a "general minimal surface which achieves a minimum 
of object dissipative energy functional. In other words, from mathematical point of 
view, this minimum problem of geometric sharp of the surface of the blade is a general 
minimal surface problem. 

Note that (4.2) is also an optimal control problem with distribute parameters, 
control variable is the surface of the blade and the Navier-Stokes equations are the 
state equations of this control problem. 

In sequel, we prove that equation (2.7) is the Eular-Lagrange equations of optimal 
problem and Navier-Stokes equations are its state equations. 

In order to investigate optimal control problem (4.2) we should consider the object 
functional J in a fixed domain in new coordinate system. 

Let WT] = ^(6)r? denote Gateaux derivative of w at 6 in direction r] e y(^2).Then 
according to (3.20) , we have following lemma 

Lemma 4.1. Assume that {w,p) is a solution of Navier-Stokes equations (3.1)(3.6) 
associated with 6 € H^{D) which define a mapping (w(6),p(6) from W{D) = 
H^{D)nH^{D) toW^iVlfxL^iVl) :Q {w{Q),p{Q)). Then {w{Q),p{Q)) is Gateaux 
differential in H^{D) with respect to all direction rj S W{D). The Gateaux derivatives 
wrj = ^^ri,pr] = ^^rj are a solution of (3.37). Then the strain rate tensor e.ij(w) 

defined by (3.2) possesses Gateaux derivatives — ^gr^(0)?? at any point @ G H^{D) 
along every direction r] G W{D), and 

^^^V = eij{w)ri + e%{w)r]x + e^^{w)r,xa, (4-3) 



Boundary Shape Control of Navier-Stokes Equations 



15 



where 



e33M=V'33H, e^fH=0. 
Proof. From (3.20) we claim 

Ral3{w, Q)t] = ^r2w'^[(5„A0fT/3 + <5/3A0a<7)»yA + (<5/3A0a + '^aA0/3)'7Aa] , 

^3a(w, e)ry = ier^w-fJ^A, i?33(e, 6) = 0. 

This leads to (4.3) and (4.4). The proof is completed. □ 

Lemma 4.2. Under the assumptions in Lemma \4-.1\ the dissipative functions <i>(w) 
defined by (4-1) is Gateaux differential at Q & H^{D) along any direction rj e W, 
and 



^-r, = $"({?, ii;)r/ + Q)r]x + $^'^(u;, e)?7Aa, (4.5) 



where 



^°{w,w) =2A'i'''e^j{w)ekiiw), 

^^{w, 9) = 2A^3>'^e^j{w)ekiiw) + Af^(w, 9), (4.6) 
<i>>"^{w, 9) = 2A'^''^e^[{w)eki{w), 

M^{w, 9) = ^e-^{Qx5a[} + <da5px)e3a{w)e3f3{w) 

+4e-4r-2aeAe33(He33(w) - Ae-^e3,a{w)eax{w) (4.7) 
-4e-3^-2(2r29«9A + (r2|V9|2 + l)(5„A)e33(u;)e3a(w^). 

Proof. Indeed, taking (3.24) and (4.3) into account, 

-^^r]e,,{w)eu{w) = M\w,e)r^x, (4.8) 
where M^{w,Q) is defined by (4.7). Then 



= iA^^f'^^^r^eM + ^^e,,(«;)e.Ku;), 
^*'"'),7 =2A'^'='e,,(iB)e„(t^)?7- 
+2A^-'"fc'e^/(u;)efc,(?«)77A, 



^77 =2A'^'='ey(zB)eH(w^)'7+(2A^^''''e<^^.(u;)+M^(?«,9))77A (4-9) 



From this it yields (4.5)(4.6). The proof is complete. □ 

Theorem 4.1. Under the assumptions in Lemma \ 3.2\ the object functional J 

defined by (4-1) has a Gateaux derivative gradj = in every direction rj ^ W :— 

H^{D) n Hq(D). Furthermore gradJ is determined by 

< gradQ{J{e)),i] >= JJ^[^"{w;w)r] + ^^{w,Q)7u + ^^''{w,e)'nx^ , ^ 
^2jir'^W'"'Qyxnxa\erdx^dx^, ^ ' ' 
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where 

$^<^(u;, e) = 2iyr^{eWf + W^O^), (4.11) 
= j^^ w'^w'^di, Wf := w"^d^. 

where w = are the Gateaux of the solution (w) of Navier-Stokes equations at 

point 6, are defined by (4-6) ■ 

Proof Indeed, taking into account of (4.5) and (4.6) we assert that 

< grade( J(e)), r? >= / j\ ^^r,erd^dx 

= J J^^ 2zy[$°(w, w)r] + $^(w, e)?7A + ^^"{w, Q)r]xa]red^dx 

D 

= J 2/i[$°(w; w)r] + ^^{w, e)r]x + ^^"{w, Q)'nxa + 2fj,r'^W'"'e^xV\a]srdx^dx'^, 

(4.12) 

where we need to show only last two terms in (4.12). In fact, from (4.4) (4.6) and 
(3.24) 

$^-(m;, 9) = 2A^^'''e^/{w)eki{w) = 4iy{e^'^^{w)eafs{w) + 2(523 ^"/S + e~^eaQ/3)e^^{w)e3(}{w) 
We^3 He33(«^) - 2e-^e0{e^''0{w)esa{w) + ea0{w)etiw)) 

+e-^eaQf3ie^'p{w)e33{w) + ea0{w)e^^{w))}, 
In view of 633 (w) = therefore, make rearrangement, 

^^^{w, e) = Ai^{e^Jj{w){ea,ij{w) - 2s-^e3e3c.H + e-^Q^Qffessiw)) 

+2e^-(u;)((.9335",a + e-2e^e^)e3^(^) - s-'epea^iw) - e-'g^'&a,e33{w))}. 

Substituting (4.4) into above leads to 

$^<^(u;,e) = 2i/r2w'^[(2e-2|ve|2 - 2g33)e^e33(w) 

+2e{g^^6x0 + s-^QxQis - e-^Qisex + \^e\''S,3x)e3f3H)] 

= 4:i^e~'^w'^[ee3x{w) - &xe33{'uj)] (4 13) 

= ir^e-^w^i^i^ + e^r^^) + \e^r\^ - ^6x.)B. 

_ i£r2^e,eA + ^e^r^w-'ex.], 

By virtue of (3.20-22) and the index (A, a) in ^^''ijxa being symmetry, 

/w'^d^w^d^ = 1 ^d^{w'^w'^)d^ = ^w^w'^l^^i^-i = 0, (by boundary conditions) 
-1 J-i 2 2 



we have 



$^-(w;, e)dC = 2ur^ J^^ W'i^^ + ^^'^H + ^ur^W^Qx. 
= 2vr'^\£Wl^ + V^X'^e^ + W^^xu]- 

The proof is complete. □ 

Taking integration by part of (4.10) and considering homogenous boundary con- 
dition for T] e W(£'),it implies 



< 



grade (J(e)), r? >= j £[9Aa(2Mr^H^''"0.A+r-$^''(w,e))-aA(r$^(w,e))+r$°(«;,w)]77(ia 
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From above discussion we obtain directly Theorem 12. II The Euler-Lagrange equa- 
tion for the extremum of J is given by: 



(4.14) 

and variational formulation associated with (4.14) reads 

Finde e V(p) = {q\q G H^{D),q\-y = Oq, |f |^ = GJsuch that 
/^{(/C"'3^'^e„/3 + ^^''{w, e))r/Aa + ^^{w, e)'r]x + i°{w,w)r]}srdx, Vr; G H^iD), 

(4.15) 

where 

l^a/SXa ._ 2iyr^W°"'Sf3X. (4.16) 

5. Second Model. Let us choice the rate of work done by the impeller under 
certain ratio of inlet and outlet pressure as minimizing functional (3.10) 



^(^) ~ / / (T ■ n ■ egrudS, (5-1) 

where n is unite normal vector to the surface 9 and {er,e0,k) are bases vector in 
rotating cylinder coordinate system, a is total stress tensor for compressible flow. 
Our purpose is to find a surface 3 of blade such that 

/(3) = if^HS), (5.2) 

where denotes a set of regular surface spanning on a given Jordannian curve C G R^. 
Under new coordinate system (5.1) can be rewritten by 

Next we compute (5.3). To do this, we give the foUowings without proof: 



(5.4) 



n{0^nHOe^, n\0^-rOx/V^, n3(-l) = (er)-i 
n3(l) = -(er)-ii+^, eg = {rey^es, (e^)" = 0, {egf = {re)-\ 

Thanks to (2.5) and (5.4) we claim that 

{-p + Xdivw)gijn'{eey = {-p + Adivw)(re)"i (53371^ + gaxn^) 

= (-p + Xdiyw){ren^ + re„n"), for ^ = ±1, (5.5) 
2^eij(w)n''{eey = 2/i(r£)-i(e33(w)n3 +e3A(w)n^), 

Substituting (5.4) (5.5) into (5.3) and simple calculations lead to 
J /(3) = /^{(l + r2e2)(-2P + 2AdivT^^) - r'^\V@\'^{-2P + 2AdivM^) 



I +2^l{re)-^{{re)-^{l + r'^QDes^iW) - rQo.e^o.{W))}rLoerdx, 
where we use notation 



(5.6) 



= w(-l), w\^^i=w{l), W ^{w{l)+w{-l))l2, W ^{w{-l)-w{l))/2, 
P = (p(l) +p(-l))/2, P = (p(-l) -p(l))/2, 

W . VW n . VP Tjcr . VW 5 . VP 

We .- Pe We .- Pe .- ^g, 

(5.7) 
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Hence 

P - P = -pi-l),W -W ^ -w{-l), (5.8) 
and in view of Lemma [3.21 the gradient of / along any direction 77 € Hq{D) is given 

by 

Theorem 5.1. Under the assumptions in Lemma \ 3.S\ the object functional I 

DI 

defined by (5.1) has a Gateaux derivative gradi = in every direction rj G Hq{D). 
Furthermore gradI is determined by 

< gradlC^),r] >= J^{[i'o{w,p,e)i] + ■i'^{w,p,e)T]x + {w,p,e)i]xa]rujer}dx, 

(5.9) 

where 

*o(^«,P,6) := (1 + r^Ql){-2P0 + 2XdivWe) - r^\VQ\'^{~2P0 + 2XdivWe) 

+2fi{re)-\{re)-\l + r^Q^essiWe) - re„e3a(W^e)), 
*^(u;,p, 8) = ■9^iw,p) + *^(u;,p)e,. + *^^(w,p)e^e^ + 'i'iiw,p, 8), 
*-^'"(w,p, 8) := -2fie-^eae^^{W) = -fir^WSx, 

(5.10) 



^^{w,p) ^2^i{{rer^i,l,{W)-s-^ip3xiW)), ^ _ 
*;)(«;, p) = [-2r^{-2P + 2XdivW)6x^ + {2r^{-P + 2XdivW) 

+26-^33 (W^)>A^2. - 2,,e-\^l{W) + r^xiW))], (5 

[ ■f^{w,p,e) =-fir^W^Qxa, 

Proof. Firstly, by virtue of (3.21) we assert 

PdivW^ , VW 

In addition, thanks to (4.7) and denote 

VW ~ VW VP ~ VP 

^^'■^VQ' ^^^=W ^«^=pe' ^^^^pe' 

(5.6) shows 

< grade/, 77 >= /o{2r202?72(-2P + 2AdivW?) + (1 + r282)(-2Fe + 2AdivW?e)?7 

-2r28A77A(-2P + 2AdivT/F) -j^\'^e\^{-2Pe + 2AdivM^e)»y_ _ 
+2^lirer'[{rer^2r^e2V2e33{W) + {re)-\l + r^Ql)^{W)r, - rrj^e^^iW) - rQ^^iW)r,] 

In view of (4.7) (4.8) and rewriting integrated we claim the first of (5.10). Moreover 

*^(w,p,8) =2r2e252A(-^ + 2Adiviy)-2r^A(-2F + 2Adiviy) ^ 

+2Ai(r£)-i[(r£)-Hl + '^'02)6^3(1^) + 2r\re)-H2x02e33{W) 
^rOo^eUW) - re^xiW)], 

(5.12) 
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Thanks to 

and Takmg (4.8) and (3.16)(3.17) into account, simple calculation from (5.12) shows 

+ [(2r2(-F + 2AdivW^) + 4fie-^^33{W))62xS2. 
^2r\-2P + 2AdivM^)<5A. - 2^JLe-\^l{W) + V^xi^mu 

+ {2^LS-\^l^{W)62. + 2^-3^3 (H^)<52a)<52^ ^ ^ 

^2iie-\r3lm + i,ll[W) + V'aA W)]©-©/^ 
-'2[i{erY^r\erW''Qxa, 

This leads to second of (5.10). 

Integrating by parts in (5.9) yields 

< Gradel,ri >^ /d[£?'^^*o(w,p, 9) 

-^{£r^uj^^{w,p,e)) + ^-|i-^(ewr2^'^-(u;,p,e))]77dx 



(5.14) 



Hence the stationary point of minimum problem should satisfies following Eular- 
Lagrange equation 

er^LO'^oiw^P, e) - £^{sr^Lu-^\w,p, 9)) + ^^{ecur^-qj^- {w,p, 9)) = (5.15) 
However, 

_|i_(^^^2^A.(^^p^0)) = -^-^{^.ujer^W'^Qx) 

= -a^i^^^er^Wnex - f,eu;r^W'^^{Qx) - ^i^^u;er^Wn^iex) 

In addition 

-£x{'fHw,P,Q)) = ^^ieiotir^WQx) = £>:{ecufir*WneaX+f^eujr^W-e^xx, 
Therefore 

-^{^i{w,p, 9)) + g-^{eu;r^^if^-{w,p, 9)) = -^^{pu^er^wnOx ~ ^{^iiver^W'')Ae, 
where 

A9 = 9ii+922. 

On the other hand, 

~dxi^^ + ^-^9, + ^';:,9,9^) = -(^A^i^ + dx'^te. + aA*^^9,9^) + (*^ + (*^^ + ^pe^)Qx. 
To sum up and by simply calculation, we obtain Theorem 12.21 

-iKoiw)Ae + K^'^iw, 9)9,a) + F"^^(w)9,9^ + F^(«;)9a + Fo{w, 9) = 0, 
9|^ = 9o, l = dD, 

(5.16) 
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Ko(w) = fxeu; ' 

K^'^iw) = a;er2(v|/^(«;) + (*^^(«;) + ^^^^{w))e^), 

f = -toerirdx-^i^iw) + 252A*^^ (5-17) 
F\w) = -ifxeudl^ir^W) + uerirdx^^iiw) + 262x'^i{w)), 
[ Fo{w, 6) = ujerHoiw, 6) - (jsr{rdx'i>^{w) + 262x'^oM), 

where ^'o,^'^ and "if^^ are defined by (5.10) and (5.11) respectively. Variational 
formulation associated with (5.16) is given by 

f Finde G V = {q\q e H^{D),q\^ = Qq} such that ,^ 

where '^o{w,p,Q),^^{w,p,Q),'<b^''{w,p,G) are defined by (5.10). □ 

6. Existence of Solution to the Optimal Control Problem. In this section, 
we discuss existence of the optimal control problem (3.1) for incompressible case. As 
well known the object functional 

j(e) = ^ J^^A'^''\e)eij{wie))eki{wie))^dxdc = ^a{w{e),w{e)), (6.1) 

is depending upon the existence of solution w to Navier-Stokes equations. Since the 
solution w of Navier-Stokes equation and are the functions of 6 , therefore J is 
a function of 8. However J can be read as a function of w : J{Q) = J(w{Q)). As 
well known that if there exists Gateaux derivative of J(0) with respect to Q at 
9*, Then the minimum point © of (3.1) is necessary to satisfies 

gradsJ{Q) = 0. 

and from this, if is regular then Eular-Lagrange equation for 6 can be obtained. 

At first, It is well known that following theorem give a sufScient condition of the 
existence 

Theorem 6.1. (Generalized Weierstrass Theorem) [13] Let X be a reflexive Ba- 
nach Space, and U a bounded and weakly closed subset of X. If the functional J is 
weakly lower semi- continuous on U,then J is bounded from below and J achieves its 
infirm on U. We consider functional defined in a closed convex set in a Sobolev space 

V{Q) := {u\u e H^'P{Q),u\ra =0,dn = To U ri,meos(ro) ^ 0}. 

Let 

J{w) = / A''''^eij{w)eki{w)y/gdxd^ 



n 

Lemma 6.1. J{w) is weak lower semi-continuous with respect to w in H^{Q)^. 

Proof. We firstly need to establish the uniform positive definitencss of the three- 
dimensional tensor A^^'''- . "Uniform" means with respect to all points x £ Q and 
to symmetric matrices of the order three {tij}. Namely, there exist a constant 
c(f2, 0, yu) > such that 

A'^''\x)tkitij>cJ2\tij\^- (6-2) 

i,3 
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See the proof of Thl.8-1 in [P.G.Ciarlet 2000]. Note that is a Lipschitz domain,^ = 
det{gij) = £^r^ > 0. Therefore, we claim 

a{w,w) > 1^0(9,8)^ ||ejj(w)||g_n 

On the other hand, the Korn inequality with boundary conditions in curvilinear co- 
ordinates (see Thl.7-4 in [P.G.Ciarlet 2000]) shows that 

Ell^^^-Wllo,^^^^("'©)IHIlo, Vi;GK (6.3) 

Hence 

a{w,w)>^lcin,e)\\w\\ln, (6.4) 

We assert that a(-,-) is a equivalent norm in H^{rt)^. In other words , J(-) is a 
equivalent norm in H^{fl)^. As well known that the norm in a Hilbert space , as 
a functional, is weakly lower semi-continuous, we assert J is weakly lower semi- 
continuous with respect to w. Proof is complete. □ 
By virtue of Lemma 13.31 we directly to obtain 

Lemma 6.2. If the function w{Q) of Navier-Stokes equations satisfies the follow- 
ing: 

Assumption P : Q„ Qoiweakly) ^ Wn — w{Qn) wo = w{Qa) (weakly). 

Then functional J(0) is weak lower semi- continuous with respect to Q. 
Finally, we have 

Theorem 6.2. Assume that the solutions {w,p) of Navier-Stokes equations with 
mixed boundary condition are weakly continuous with respect to Q.Then there exists 
a two dimensional surface S defined by a smooth mapping 

e-.D — >w = H^{D) n (6 e h^{d),q\qd = e*, ) 

such that J(0) achieves its infirmum at {Q.,w{Q)}. 

Next we consider the existence of the solutions for Navier-Stokes equations. In 
deed flow's domain is unbounded domain. In section 3 we make artificial bound- 
ary, inflow boundary Tin and outflow boundary Tout, and impose natural boundary 
conditions (3.6). We also can impose the pressures 

P\ri„ ^Pin, p|r„„t =Pout, 

or flux 

/ pw ■ ndT = Q, / pw ■ ndT = Q, 

Let us consider energy inequality. Owing to 

(2uj X w,w) — 0, 
Hence moment equations (3.11) show 



a{w, w) -\- b{w, w, w) — (/, w), 
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However, 

b{w,'w,'w) = j^w^V jw"^ gikW^ ^/gdxdS^ = J^{Vj{w^w^) — w^divw)gikw'^y^dxd£, 
= J^{diy{\w\'^w) — QikW^w^W jw'^)^dxd^ = Jp^ \w\'^w ■ ndT — b{w, w, w) 



b{w,w,w) = - / \w\'^'w-ndr, (6.5) 
2 Jti 

Here we denote 

\wf = Qikw'w^, Ti = U Tout- (6.6) 
The flux of inflow and outflow kinetic energy are respectively give by 

Kin{w) = / \w\^w ■ ndT, Kout{w) = / \w\'^w ■ ndT 

where w ■ n = gijW^n^ and n is outward normal unite vector to inflow or out flow 
boundary. Therefore (6.1) shows that 

h{w, w, w) = Kout{w) - Kin{w), (6.7) 

It is obvious that we can not make say Kout{w) < oo . 

Theorem 6.3. Suppose that the exterior force f and normal stress g of inflow 
and out flow boundaries Fi = Tin U ^out satisfy 

\\F\U ll./llo,o + ||.9m||-i/2,r,„ + ||,go«dl-i/2,r„„t < i^V(4ciC2), (6.8) 
and the mapping Q defined in (2.1) is C^{D)- function satisfied by 

Ge C^D,Ko) = {eG <^^.^||^^^|| <^^|^ (g_9^ 

where Kq is a constant. Then there exists a smooth solution of the variational problem 
(3.11) satisfying 

where ci,C2 are constants depending of{Q,Kct) defined, by (6.14)(6.15). 

Proof. To prove the theorem for a steady Navier-Stokes problem, it is conve- 
nient to construct the solution as a limit of Galerkin approximations in terms of the 
eigenfunctions of the corresponding stead Stokes problem. This use of the Stokes 
eigenfunctions (see in Heywood and Rannacher and Turek[20] and Glowinski[12]). 
Galerkin equations are a system of algebraic equations for constant unknowns and 
the Galerkin approximation solution w is a solution of the flnite dimensional problem 

a{w, v) + 2{ii)xw, v) + b{w, w, v) =< F,v >,yv G V"to := span{<^i, (^2 • • • , <^m} (6.11) 

where (pi,i = 1,2,- ■ - m are eigenfunctions of corresponding Stokes operator. Let Sp 
denote the spheres in Vm satisfying inequality (6.10). Assume that G Sp, we 
seeking w such that 

a{w, v) +2{u> X w,v) +b{w^,w,v) =< F,v >,\/v & Vm, (6.12) 
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(6.12) is uniquely solvable because S Sp because the w = is only solution of the 
corresponding honiogcncous equation (F = 0). Indeed, if satisfies (6.10) and w 
satisfies (6.12) with the F = 0, taking (6.4) and 

2{(jxw,w)=0 (6.13) 

into account, and combing (3.19) and (6.11) we assert 

c{fl,Q)n\\Vw\\l^< \b{w^,w,w)\ < ci(f2,/<:o)||w,||i,6||Vw||o||w||2 < ci||Vu;*||o||Vw||g 
<ci(0,i^o)^^||VHI^, 

(6.14) 

This implies that 10 = 0. In order to apply Brouwer's fixed point theorem, we have 

to show the mapping w; take the ball Sp defined by (6.10) into itself, suppose 

that satisfies (6.10). Similarly to (6.13), for nonhomogeneous equation {F ^ 0) 
,we obtain 

c{n,e)n\\Vw\\la < \b{w,,w,w)\ + \<F,w>\< 
ci{Q,Ko)\\Ww4o\\VwU + C2iQ,Ko)\\F\\4Vw\\o, 

c{n,e)fi\\Vw\\o,n < ci(0,i^o)||V«;,||o||V«;||o + C2(0,i^o)||i^||*, (6.15) 
Here note that 

Fw = gijF'uP = {6a0 + r'^eaQ0)F°'w'^ + eQaiF^w^ + F^w") + eV^F^w^ 
Therefore 

llyy^ll < C2\\F\\, ^ C2\\F\\, 

" " - c(o,e)M- ci||v«>||n - c(n,e)M-^^^^ ^^[i-^(i-4ciC2||F||./(c(a,e)j:^)^)] 

Thus Brouwer's fixed point theorem can be applied and gives the existence of Galerkin 
approximations satisfying (6.8). Hence by a standard compactness argument there is 
at least a subsequence of the Galerkin approximation converging to a weak solution 
w e V{Q,) of the steady problem (3.11): 

Fmd{w,p),w e V{n),pe such that 

a{'w, v) + 2{uj X w,v) + b{w, w, v)+ 

-{p,divv) =< F,v >, Vw e V{n), 
(g,divw) =0, Vg e i.=^(0). 

Its smoothness is easily proven if one obtain a further estimate from the Galerkin 
approximations by setting v = — Aw in (6.9). This gives 

/imwllo = -2{(j X w, Aw) - b{w, w, Aw)+ < F, Aw >, (6-14) 

Because Aw is solenoidal, one has the rather unusual trace estimate 

\2{u;xw,Aw)\<C3\\wU\Aw\\n, \ < F, Aw > \ < C3\\F\\4Aw\\o, (6.15) 
which we combine with (6.14) and Agmon inequality 

Iklloo < c4Vw\\l^^\\Aw\\l/\ yw e DiA), (6.16) 
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to get 

^J-\\Aw\\l < c4Vw\\l^^\\Aw\\l^^ + C3\\w\\o\\Aw\\o + C3||F||,|| A«;||o. (6.17) 
Then, by using Young's inequality, we obtain 

f^WAwWo < ^-§\\^w\\l + ^IIHI^ + ^ll^ll^ (6.18) 

which is then inherited by the solution. The full classical smoothness of the solution 
can now be obtained using the i^-regularity theory for the steady Stokes equations. 
This completes the proof of Theorem 16.31 □ 

It is clear that the bound in (6.10) is not uniform with respect to 8. Next our 
gaol is to prove the the solution w(8) of Navier-Stokes equations (3.11) with mixed 
boundary condition is uniformly bounded. At first we prove 

Lemma 6.3. Assume that 

Be T={<Pe C\n),sup |V</)| < K, := -^,} (6.19) 

D Zy'br 

then 

a{w{Q), w{&)) > n/2{2eaf3iw)eaf3{w) + e3a(w)e3a(w) + e33(w)e33(w)), Vw e V{n) 

(6.20) 

where a(-, •) is defined by (6.1) and \\ ■ \\n is defined by (3.25). 
Proof. By virtue of (4.6)(2.1)(A — 0) we claim 

a{w, w) = 2ii[eai3(w)eai3{w) + {2e-'^QxQa + 26xag^^)e3\{w)e3a(w) + g^V^e33(u;)e33(u;) 
-A£-^Qae3\{w)ea\{w) + 2£-'^QxQae\a{w)e:i:i(w) - 'ie~^Qxe?.x{w)e:i3{w)], 

(6.21) 

By Cauchy inequality and Yang inequality we obtain 

4£-ieAe3A(t«)e33(w) < 8e-2(e^e3A(H)' + W^^^(w)f^ 



4£-le„e3A(w)eaA(w) <4 /£-2^(e„e3A(u'))2 /E(eaA(w^))2 

< 16£-2|ve|2(e3A(w)e3A(if^)) + je^x{w)eo,x{w), 
2£-^QxQaexa{w)e3z{w) < 2^£-^j:iexe.e33{w))^^j:iexAw)y 

< 4£-4|ve|4e33(w)e33(ti') + ieAa(w)e„A(w), 

(6.22) 

To sum up and thanks g^^ = e^^r^^ + £^^|V8p we assert that 



a{w, w) > 2n[^ec0{w)ec,p{w) + [-6e-28A8^ + 2e-^r~^Sxa ~ 14e-2| v8|2,5A^]e3A(«^)e3<T(w) 
Note 



+ (i5'V' - 4£-4|V8|4)e33(He33(w;)], 



8Ae3A(^i;)e.e3.(^i') < (JEQI/E^3aM)' < |V8|2e3A(«;)e3A(w) 



(1^33^33 _ 4,-4|ve|4) ^ + 2r-2|V8p - 7|Vn 
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Therefore 

aiw,w) > 2^[|e„/5(w)e„/3(w) + 2e-'^r-^[l - 10r2|Ve|2]e3A(w)e3AM 
+ ^^(1 + 2r2|Ve|2 - 7r4|Ve|4)e33(z«)e33H] 

It is obvious that 

2£-V-2[i _ 10r2|Vep] > if |Ve| < 

^^^(1 + 2r2|Ve|2 - 7r4|Ve|4) > \e-^T-^, if |Ve| < ^, 

Hence If 9 G then 

a{w,w) > 2fi[^eaf3iw)eaf3iw) + je3x{w)e3x{w) + 1633(^)633 (w)], 

End the proof. □ 

Lemma 6.4. Assumption in Lemma \4-.l\ is held such that the constant Kq satisties 

Ko < , (6.23) 

where Poi^) 'is a constants depending Q only .Then we have following unform coer- 
civeness for bilinear form a(-, •) 

a{w,w) > ^\\w\\l > i^C^m\w\\l,,, Vw e H\nf, (6.24) 

where ao{i^) is a constant independent of Q and is defined by (3.25). 

Proof. In view of (3.16), 

\Mw)\\l,, = \Wc.p{w)\\ln + ll^;^Mt^)eA||g.o + U^^%w)QxQ4l^ + \Kp{w)\\l^ 

+2{^^p{w), iPlp{w)Qx) + 2((p„0(u.), V^;^(u;)eAe,) 
+2[^l%{w)Q,Q^Ji, ^^^p{w)Qx) + 2(e„;3(^), e:^(u;)), 

(6.25) 

Using (3.18) it is easy to verify 

Ulp{y^)\\ln < Poim^MWl^, \K''p{w)\\l^ < (3om{y^^iw)\\i,, + ||^33(HII§.o), 

uLMWln < MmyMw)\\in + y33{w)\\in) 

ll</3(^)llo,n < PomK'o\Mw)\\o,n, \Mw),el^{w))\ < K^/^poimwlH, 

(6.26) 

Remark 6.1. the constant Pq{Q) represents different meaning at different place. 
By means of Schwarcs inequahty, from (6.25) (6.26) we claim 

\Kf3{w)\\l^ > \Mw)\\l^ ~ Kl^'PomMh, (6.27) 
By similar manner we imply 

||e3o(Hf > y3c.iw)\\l,,-K'J'Pom\Ml, WessHW^ > \Mw)\\l,,~KoPomMl, 
Let return to (6.20), it yields that 

aiw,w) > >l[\\w\\l,, - K'J^MmMln]- 
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This yields that if (6.23) is satisfied then ,by Lemma[331 it yields (6.24) . □ 
From Theorem 15.11 and Lemma 14.21 we conclude 

Theorem 6.4. Assume the the assumptions in Theorem \ 6.1\ are satisfied and 
(6.23) held. The w{Q) is a solution of rotating Navier-Stokes equations (3.11) asso- 
ciated with Q E C'^{ fl, Kq). Then the following estimate held 

Me)||o<|^[l-^(l-^^2ig5|^)], (6.28) 

The solution is sequently weak continuous with respect to Q, i.e. if sequence {Qi} 
is weak convergent in H^{D) then there exist convergent subsequent in {w{Qi)} in 
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